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IVo  possible  sources  of  error  in  computing  a tropospheric  range  correction  from  a model  are 
a mismatch  between  the  model  and  the  actual  refractlvity  profile,  and  neglect  of  a signal  path  bend-1 
Ing.  To  study  the  error  sources,  the  tropospheric  range  effect  on  a radio  signal  has  been  evaluatedj 
at  different  elevation  angles,  for  some  observed  atmospheres,  by  three  methods:  (a)  using  meteoro- 

logical balloon  data,  the  effect  is  evaluated  along  the  computed,  curved  signal  path  (this  is  our 
best  estimate  of  the  true  effect  and  is  the  standard  for  comparison);  (b)  using  meteorological  bal- 
loon data,  the  effect  is  evaluated  for  assumed  straight-line  propagation;  and  (c)  using  only  surface| 
meteorological  data,  the  effect  is  evaluated  from  the  two-quartlc  tropospheric  model  for  assumed 
stralghc-line  propagation.  The  differences  ((b)  - (a))  and  ((c)  - (b)),  respectively,  represent  the! 
path-bending  component  and  the  profile-mismatch  component  of  the  total  model  error  ((c)  - (a)).  Re-| 
eults  are  presented  in  figures.  The  effect  of  the  model-based  correction  on  the  computation  of 
tracking-station  poaltiur.  is  shown  fot  a few  sample  cases.  ^ 
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1.  INTRODUCTION  AND  BACKGROUND 


Tropospheric  effects  on  satellite  range  data,  and  the  diffi- 
culty In  predicting  them,  are  greater  at  low  elevation  angles  than 
at  high  because  of  the  Increased  length  of  the  signal  path  in  air, 
especially  In  the  lower  layers  of  air.  The  apparent  electromag- 
netic range  to  a satellite  Is  jnds  along  the  signal  path,  where  s 
Is  the  signal  path  length  and  n is  the  Index  of  refraction  of  the 
medium  and  varies  along  the  path.  The  atmospheric  effect  on  the 
measurement  Is 


curve 


Jnds  - p , 


where  p Is  the  slant  range.  If  there  Is  no  significant  curvature 
of  Che  path,  the  integral  can  be  taken  along  the  slant  range  vector 
and  the  effect  can  be  written 


sc.  line 


- J(n  - l)dp  . 


For  convenience,  the  refractlvlty  (N)  Is  defined  as 

N - 10®(n  - 1)  . (3) 

A detailed  knowledge  of  the  refractlvlty  profile  at  the  per- 
tinent time  or  a suitable  model  of  It  Is  necessary  to  evaluate 
Eqs.  (1)  or  (2)  and  so  to  correct  for  the  tropospheric  effect. 
Detailed  knowledge  of  the  profile  is  seldom  available.  Any  reason- 
able model  can  be  adjusted  to  give  the  proper  Integral  of  Eq.  (2) 
at  the  zenith,  but  only  a model  that  also  gives  the  correct  pro- 
file shape  will  give  the  correct  integral  of  Eq.  (2)  at  both  high 
and  low  angles.  The  effect  predicted  by  the  two-quartic  refrac- 
tlvlty profile  model,  will  be  compared  with  and 

^st.  line* 

Path  bending  and  model  inaccuracy  have  greater  effects  near 
the  horizon  than  at  high  angles.  This  study  examines  both  of  these 
lov-angle  effects. 

Refractlvlty  is  a function  of  meteorological  conditions. 
Smith  and  Ueintraub  (Ref.  1)  give  the  following  expression,  for 
radio  frequencies  up  to  at  least  15  GHz; 

Ref.  1.  E.  K.  Smith,  Jr.,  and  S.  Welntraub,  "The  Constants 
In  the  Equation  fov  Atmospheric  Refractive  Index  at  Radio  Frequen- 
clos."  Proc.  I.R.E..  Vol.  41,  No.  8,  August  1953,  pp.  1035-1057. 
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(p  ^ ) . 


(4) 


Here  T is  the  Kelvin  temperature,  and  P is  the  total  pressure 
and  e the  partial  pressure  of  water  vapor,  both  in  millibars. 

The  above  two  terms  of  N are  the  "dry"  and  "wet"  terms,  respec- 
tively (subscripts  d and  w).  The  two-quartic  refractivity  pro- 
file expression  was  developed  from  a study  of  meteorological  bal- 
loon data.  It  deals  separately  with  the  height  variation  of  the 
tv;o  terms.  The  full  expression  is  (Ref.  2): 


N 


- + N 

d V 


N 


(h  - h )**  ^ 

d 6 


(h^  - h)**  h £.  h. 


(5) 


(h  - h ) 
w s 


(h  - h)‘ 


4 W 


h S h 


In  these  equations,  h is  height  above  the  geold  and  the  subscript 
s refers  to  surface  values.  The  heights  hj  and  have  been  de- 
termined from  balloon  meteorological  data  (Refs.  3 and  4). 

An  expression  for  Apmojei  at  an  elevation  angle  E has  the 

form 


Ap  , , = 10  ^ 
model  ^ 

i 


G(E,h.) 
«i  ^ 


(6) 


Ref . 2 . H.  S.  Hopfleld,  *'Two-Quartic  Tropospheric  Refrac- 
■tlvlty  Profile  for  Correcting  Satellite  Data,"  J.  Geophys.  Res., 
Vol.  74,  No.  18,  20  August  1969,  pp.  4487-4A99. 

Ref.  3.  H.  S.  Hopfleld,  "Tropospheric  Effect  on  Electro- 
magnetlcally  Measured  Ranges  Prediction  from  Surface  Weather 
Data,"  Radio  Sci. . Vol.  6,  No.  3,  March  1971,  pp.  357-367. 

Ref.  4.  H.  S.  Hopfleld,  "Tropospheric  Range  Error  Parame- 
ters: Further  Studies,"  APL/JHU  CP  015,  June  1972  (also  published 

as  Goddard  Space  Flight  Center  Preprint  X-551-72-285,  August  1972). 
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where  1 takes  the  values  1 and  2 for  the  dry  and  wet  components, 
respectively.  The  full  expression  for  the  function  G(E,hi)  is 
given  in  Ref.  2.  A simplified  expression  for  easier  computing  is 
given  in  Ref.  5. 

The  following  assumptions  and  approximations  entered  into 
the  development  of  the  model  and  the  correction: 

1.  Ray-tracing  methods  can  be  used  to  determine  Che  signal 
path;  Interference  and  diffraction  phenomena  are  not 
considered.  This  assumption  is  retained  here. 

2.  The  temperature  lapse  rate  in  the  un-ionlzed  atmosphere 
is  approximately  6.8®C/km  and  is  constant.  This  is  a 
good  approximation  for  the  troposphere  proper  but  not 
for  the  stratosphere;  however,  the  h'd  model  lumps  them 
together.  Since  807,  of  their  joint  effect  occurs  in 
Che  troposphere  proper,  below  the  tropopause  (a  still 
larger  percentage  at  very  low  angles) , the  resulting 
mean  error  in  using  the  model  may  still  be  minor. 
Deviations  from  the  mean  may  be  significant,  however, 
since  temperature  inversion  layers  and  ocher  irregu- 
larities actually  are  frequent  but  are  disregarded  in 
the  model.  Their  effects  will  be  seen  in  this  study  but 
only  on  a sample  basis. 

3.  The  N„  profile  is  written  in  the  same  quartic  form  as 
the  Nj  profile  but  does  not  have  the  same  theoretical 
justification.  Observed  profiles  are  very  Irregular. 

A.  Curvature  of  the  signal  path  is  neglected  in  deriving 
the  correction  based  on  the  model.  This  effect  is  ex- 
amined here. 

5.  Horizontal  gradients  of  refractivity  are  still  ne- 
glected in  the  present  study;  weather  fronts  are  not 
considered. 

The  study  is  aimed  toward  improving  the  tropospheric  cor- 
rection for  low-angle  range  and  range-rate  data.  A little  addi- 
tional theory  has  been  developed.  Items  2,  3,  and  4 in  the  above 
list  of  approximations  have  been  studied  with  the  help  of  some 
sample  data.  Various  interesting  findings  are  presented  that  have 
not  yet  been  assimilated  into  a final,  unified  correction  for  low- 
angle  data.  This  document  is  a progress  report,  not  a final  re- 
port on  the  Subject. 


Ref.  5.  S.  M.  Yionoulis,  "Algorithm  to  Compute  Tropospheric 
Refraction  Effects  on  Ranee  Measurements,"  J.  Geophys.  Res..  Vol. 
75,  No.  36,  20  December  1970,  pp.  7636-7637. 
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2.  THEORETICAL  CONSIDERATIONS 


PROCEDURES  FOR  STUDYING  LOW-ANGLE  EFFECTS 
Ray  Tracing 

The  comparison  of  AP curve  line  depends  on  knowing 

the  path  taken  by  a signal  traveling  through  the  ac-nosphere. 

In  collaboration  with  Harry  K.  Utterback,  a computer  pro- 
gram was  developed  for  tracing  the  cui'^ing  ray  path  cf  a signal  at 
any  elevation  angle  through  an  atmosphere  with  any  specified  re- 
fractivity  profile  (see  Fig.  1).  Details  will  be  published  in  a 
document  that  is  now  being  prepared.  In  brief,  the  program  com- 
putes the  geometrical  length  of  the  curved  path  and  also  its  elec- 
tromagnetic length  (Jnds),  using  numerical  integration.  Both  the 
geometrical  and  the  electromagnetic  path  lengths  are  compared  with 
the  slant  range  vector  (p)  having  the  same  endpoints.  The  total 
bending  angle  is  also  determined.  The  program  can  use  either  a 
theoretical  refractlvlty  piofile  or  one  computed  from  observed 
meteorological  balloon  data  (Ref.  3).  Observed  profiles  are  used 
here. 

The  hc-*2hts  at  wM<'h  meteorological  data  are  measured  by 
balloon  Instruments  are  considered  as  boundaries  of  successive 
atmospheric  layers.  The  zenith  angle  (Bj)  of  the  ray  at  the  sur- 
face of  the  earth  is  specified.  The  zenith  angles  (Bi)  at  all  the 
layer  boundaries  are  then  obtained  from  Snell’s  law  for  spherical 
surfaces: 

n^Cj^  sin  6^^  • constant  ("’) 

where  values  of  r^^  are  the  radii  of  the  layer  boundaries  from  the 
center  of  the  earth,  and  the  constant  is  obtained  from  the  known 
values  at  the  surface.  The  ray  path  in  each  layer  is  assumed  to 
be  an  arc  of  a circle.  There  are  no  discontinuities  In  the  ray 
direction,  but  the  radius  of  curvature  changes  at  the  layer  bound- 
aries. Each  radius  of  curvature  can  be  computed  from  the  known 
zenith  angles  and  layer  heights.  A mean  value  of  n for  each 
layer  Is  used  in  the  numerical  evaluation  of  Jnds  along  the  curved 
path. 

It  should  bfc  noted  that  the  boundary  r^^p  (Fig.  1)  repre- 
sents the  topmost  height  of  the  specific  set  of  balloon  observa- 
tions, not  the  top  of  the  atmosphere.  Special  methods  developed  to 
estimate  the  increments  of  range  effect  and  signal  path  bending 
above  r^op  will  be  described  in  the  forthcoming  paper  on  the  ray- 
tracing  program. 


10  - 


THE  JOHNS  HO*>K-N5  JNtVCRS'lV 

APPUEO  PHYSICS  LABORATORY 

LAUBEC.  MABtlASO 


The  resulting  value  of  for  a particular  atmosphere 

and  elevation  angle  is  our  best  estimate  of  the  actual  tropospheric 
range  error  under  the  specified  conditions. 

In  addition,  the  refractlvity  is  integrated  along  the  slant 
range  vector  in  two  ways:  (1)  using  the  same  refractlvity  profile 

as  the  ray  trace  to  get  APgt.  line*  ^2)  using  the  two-quartic 

refractlvity  profile  based  on  surface  data  only  to  get  ApmQjjei. 

Effect  of  Tropospheric  Correction  on  Ground  Station  Position 

A tropospheric  correction  based  on  the  above  theoretical  re- 
fractlvlty  profile  is  routinely  applied  to  doppler  data  from  navi- 
gational satellites.  The  effect  of  this  correction  on  "navigating" 
a station  position  was  studied  for  a few  sample  cases.  The  station 
position  was  determined  with  and  without  the  use  of  the  correction. 
The  effect  of  changing  the  value  of  a height  parameter  of  the  model 
was  also  examined  for  a small  sample  of  satellite  passes. 


SIGNAL  PATH  BENDING 

For  correcting  range  and  doppler  data,  the  angular  bending 
of  the  signal  path  is  of  concern  only  insofar  as  it  affects  the 
electromagnetic  path  length.  The  actual  signal  path  of  length  s 
is  usually  curved  except  in  the  zenith  direction  and  is  geometri- 
cally longer  than  the  slant  range,  but  it  is  electromagnetically 
shorter  according  to  Fermat's  principle  (Ref.  6): 

s > p I 

but  i (8) 

i nds  < J , , ndp  . ) 

‘’curve  *'st.  line  ' 

The  ray-tracing  results  show  that,  to  a good  approximation, 

jnds  - jndp  -(s  - p)  . (9) 

The.  quantity  (jnds  - jndp  ) , the  curvature  error  in  the  tro- 
pospheric correction,  must  be  evaluated  by  ray-tracing  techniques, 
but  the  quantity  (a  - p),  which  is  nearly  equal  to  it  in  magnitude, 
can  be  studied  geometrically  for  simple  cases. 


Ref,  6,  G,  Joos,  Theoretical  Physics,  translated  by  Ira  M. 
Freeman,  Hafner  Publishing  Co.,  New  York,  193A. 
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The  bending  effect  on  the  signal  path  length  is  a few 
millimeters  or  less  at  elevation  angles  of  20“  and  greater,  a few 
centimeters  at  10“  elevation,  and  a few  meters  at  1“  elevation. 

The  usefulness  of  low-angle  data  for  geodetic  purposes  therefore 
depends  partly  on  our  ability  to  correct  adequately  for  signal 
bending  effects. 

Amount  of  Bending 

Path  bending  (t)  is  produced  by  a refractive  index  gradient 
normal  to  the  path,  and  the  rate  of  bending  (di/ds)  along  the  path 
is  proportional  to  this  gradient.  In  a spherically  layered  atmo- 
sphere, the  bending  rate  is  (Ref.  7) 

, ^ sin  6 

dr  dr  (10) 

ds  n 


where  the  definitions  are  the  same  as  above.  The  negative  sign  is 
chosen  to  make  t positive  for  a path  that  is  concave  downward,  as 
in  the  real  troposphere.  Obviously  there  is  no  bending  at  the 
zenith,  where  = 0.  The  total  bending  through  the  atmosphere  is 


di 


dT  ds 


J dl  = J h7 


ds  dr 


From  geometry, 
ds  1 

- - i C3  ■■  ■ ^ 

dr  cos  $ 


Combining  equations. 


T = 


_ dn  sin  6 
**  n HT  cos  6 


.j- 


(11) 


(12) 


(13) 

(14) 


This  equation  is  well  known  (Refs.  8 and  9)  but  is  not  easy  to 
evaluate.  It  can  be  numerically  Integrated  using  a model  of  the 


Ref.  7.  J.  E.  Freehafer,  "Geometrical  Optics,"  Vol.  13, 

Chap.  2,  Propagation  of  Short  Radio  Waves,  D.  E.  Kerr  (Ed.), 

M.  I.  T.  Radiation  Laboratory  Series,  McGraw-Hill  Book  Co.,  1951, 
and  Dover  Publications,  New  York,  1965. 

Ref.  8.  V/.  M.  Smart,  Text-Pook  on  Spherical  Astronomy, 

Cambridge  University  Press,  Cambridge,  England,  1962. 

Ref.  9.  B.  R.  Bean  and  E.  J.  Dutton,  "Radio  Meteorology," 

National  Bureau  of  Standards  Monograph  92,  1 March  1966. 
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n profile  or  data  from  an  observed  profile;  or  expansions  and  ap- 
proximations may  be  used  to  get  a general  expression.  Such  approxi- 
mations are  generally  valid  at  high  elevation  angles  but  not  ade- 
quate at  low  angles. 

A new  approximate  procedure  is  described  below.  It  appears 
to  be  satisfactory  down  to  a lower  elevation  angle  than  is  ordi- 
narily the  case. 

For  this  new  procedure,  we  use 


sin  6 


nr 


sin  6^ 


(15) 


from  Snell's  law  (Eq.  (7));  the  subscript  1 now  refers  to  values 
at  the  surface  of  the  earth.  Then 


cos  6 “ 


sin^  6, 


(16) 


Dividing  Eq.  (15)  by  Eq.  (16)  and  simplifying  but  making  no  ap- 
proximations, we  get 


Vi 

nr 


tan  6 “ tan  6, 


Substituting  Eq.  (17)  into  Eq.  (14)  yields 


i 1 + tan^  6, 

1 1 

L V nr  / 

T “ -tan  6 


J 


Vi 

nr 


dn 


(17) 


(18) 


In  the  troposphere  and  stratosphere,  n decreases  as  r increases. 
The  product  nr  therefore  changes  more  slowly  than  either  factor 
separately.  The  ratio  n^r^/nr  changes  by  less  than  1%  up  to  50  km 
above  the  earth  and  by  less  than  0.1%  in  the  lowest  5 km  of  air, 
where  half  of  the  path  bending  occurs.  It  may  be  satisfactory  to 
assume  a constant  value  for  this  ratio,  i.e.,  some  reasonable  mean 
value. 


14 


TH6  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREt.  MARYLAND 


On  this  basis,  let 
nr 

R = ■ ■■  ■ = constant 
nr 


and  let 


^ 1 + tan^  8^  (1  - R^) 

which  is  no  longer  a function  of  r.  Using  these  values,  Eq. 
becomes  $ 

1.0 

T = - A(B^)  tan  \ J ^ * 

"l 

•Oiis  Equation  can  be  integrated  to  yield 
T = A($j^)  tan  An  n^^  . 

For  2 > n^  > 0 , 


(20) 

(18) 


(21) 


(22) 


An  nj^  = (n^^  - 1)  - i(nj^  - 1)^  + s(n^  - 1)  ^ . 


Since  (n;j^  - 1)  is  a very  small  quantity,  only  the  first  term  need 
be  used.  Substituting  a value  from  Eq.  (3),  we  have  for  the  total 
bending  angle 

T = A(B^)  tan  B^  x 10“®  rad  . (23) 

This  equation  cannot  be  used  when  6i  = 90®  (horizontally  arriving 
ray).  However,  it  does  not  imply  an  infinite  value  of  t,  since 
A(B2^)  ->■  0 as  tan  6i  -*■ 

Values  of  t computed  from  this  equation  will  be  compared  in 
a later  section  with  values  obtained  from  ray  tracing. 


Range  Effect  of  Bending 

This  effect  is  easily  deduced  for  a hypothetical  case  where  the 
path  of  the  signal  is  an  arc  of  a circle  (see  diagram) . 
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The  signal  travels  along  a circular  arc  of  length  s to  or 
from  an  object  at  a slant  range  p and  is  deflected  through  an  angle 
T.  The  radius  of  curvature  is  r^,.  From  the  geometry, 


s = r T 
c 


But 


Then 


p = 2r  sin  -r  . 
c I 


sin  1 = 1 . (^/2)  + (T/2)  _ 

2 2 3!  5! 


» • 


2r 


8 - P 


/ T_  jr 

A 2 ■ 4 


48  3840 


\ 24  1920  *'*/ 


For  small  values  of  t,  to  first  order, 


s - p 


(24) 


In  the  atmosphere,  short  segments  of  the  path  may  be  approxi- 
mated Individually  as  circular  arcs,  but  the  path  as  a whole  is  not 
circular  since  the  gradient  of  the  refractive  index  is  not  constant 
(cf.  Eq.  (10)).  No  theoretical  expression  has  been  found  to  relate 
(a  - p)  to  T for  the  real  atmosphere,  but  an  empirical  relation  will 
be  illustrated  in  a later  section. 
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3 . NUMERICAL  RESULTS 


RAY  TRACING 

Figure  2 shows  the  magnitude  of  the  tropospheric  radio-range 
effect  as  a function  of  signal  elevation  angle  at  the  surface  of 
the  earth.  Any  correction  must  try  to  match  such  effects.  The 
values  shown  are  means  obtained  from  a one-month  sample  of  meteo- 
rological balloon  data.  The  ray-traclng  procedure  was  used  with 
the  data  of  each  balloon  flight  to  determine  Jnds  along  ^the  curv- 
ing ray  path  at  the  time  of  the  flight.  The  quantity  (Jnds  - p) 
represents  our  best  estimate  of  the  tropospheric  effect  actually 
encountered  by  the  signal.  The  slant  range  in  each  case  is  the 
range  to  the  endpoint  of  that  particular  ray  trace,  1000  km  above 
the  earth.  The  integral  j (n  - l)dp  computed  along  the  slant  range 
vector  t hr  ugh  the  same  atmosphere  is  the  error  estimate  to  the 
same  endpoint  when  curvature  of  the  signal  path  is  neglected. 

Both  sets  of  results  are  shown.  The  difference  between  the  two  is 
the  error  incurred  by  neglecting  the  bending.  It  is  imperceptible 
at  high  angles  but  is  significant  at  the  horizon.  The  straight- 
line  path  gives  the  larger  tropospheric  effect  (Fermat's  princi- 
ple) . 

Figure  3 shows  the  difference  on  a larger  scale,  for  low 
elevation  angles  only.  A radio-range  correction  that  is  perfect 
at  the  zenith  but  neglects  the  signal  path  bending  is  3 cm  too 
large  at  10*  elevation  for  this  set  of  data,  19  cm  too  large  at  5®, 
and  3 m too  large  at  1®  (elevation  angle  of  arrival;  elevation  of 
the  corresponding  slant  range  vector  is  0.4®).  The  values  are 
means,  with  error  bars  at  ±o  for  the  lowest  angles.  Radio  effects 
in  summer  would  be  a little  larger;  effects  on  optical  data  (i.e., 
laser)  would  be  slightly  smaller  than  those  in  Fig.  3. 

Results  from  the  ray-tracing  procedure  provide  interesting 
sidelights  on  the  problem.  Some  of  these  will  be  included  among 
the  following  figures. 

Figure  4 shows  the  cumulative  bending  effect  on  a radio  sig- 
nal going  out  from  the  earth,  as  a function  of  height  in  the  atmo- 
sphere, for  two  sample  atmospheres,  a summer  and  a winter  sample. 
Most  of  the  bending  occurs  in  the  lower  part  of  the  atmosphere. 

J he  greater  radio  signal  bending  by  the  summer  atmosphere  is  a re- 
sult of  the  greater  water-vapor  content.  Irregularities  in  the  at- 
mosphere have  a visible  effect,  but  it  is  not  conspicuous  in  this 
small-scale  figure.  The  irregularities  are  emphasized  in  Fig.  5. 
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Fig,  3 Effect  of  Signal  Path  Bending  on  Radio  Range  Measurement , APj^  Ij^g  curve 


Height  (km) 

Fig.  4 Cumulative  Angular  Bending  of  Radio  Signal  Path  versus  Height  above  the  Earth 
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Figur*  5 shows  the  radius  of  curvature  of  the  signal  path, 
layer  by  layer  of  the  atmoephate,  for  one  of  the  ray  traces  of 
Fig.  4,  also  at  an  elevation  angle  of  10°.  The  effect  of  irregu- 
lar wacer'Vapor  and  temperature  profiles  in  the  lower  atmosphere  is 
very  conspicuous.  Above  the  water -vapor  region,  the  radius  of  cur- 
vature Increases  along  a smooth  but  not  a simple  curve  and  shows 
the  effect  of  the  change  of  temperature  regime  at  the  tropospause. 
In  one  thin  layer  near  the  earth,  the  radius  of  curvature  of  the 
signal  waa  littls  greater  than  the  radius  of  the  earth;  at  a height 
of  22  km  abova  the  earth,  it  was  aa  great  as  the  earth-moon  dis- 
tance. 

Flgura  6 ahowa  the  aame  bending  angles  as  in  Fig.  4,  but  as 
a function  of  the  rcfractivlty  along  the  path  (here  plotted  so  as 
to  daersase  to  tha  right  on  the  graph,  i.e.,  upward).  The  cumula- 
tive angular  bonding  for  a ray  at  a given  elevation  angle  is  nearly 
a linear  function  of  refractivlty.  This  is  *n  agreement  vlch  the 
theory  diacusasd  above.  The  two  atmospheric  samples  yield  paral- 
lel linos  so  chut  the  larger  total  bending  results  from  the  larger 
surface  vsluo  of  N (cf.  Eq.  (23)).  This  statement  is  illustrated 
by  the  next  figure. 

Figure  7 shows  the  total  angular  bending  for  winter  and  sum- 
mer ssmplas,  but  for  s lower  angle  than  in  Fig.  6,  plotted  against 
^surfacs’  optlcaj  (N^ry  ®loTie)  and  radio  (total  N)  effects 

ars  included  in  the  figure,  with  the  bending  angle  obtained  from 
ray  tracing.  The  data  separate  themselves  into  natural  groups, 
but  all  lie  along  a nearly  straight  line,  as  predicted  by  Eq.  (23). 
Optical  signals  (dry  component)  are  bent  less  in  summer  chan  in 
winter  (lower  air  dsnalty  gradient),  but  radio  signals  are  bent 
mors  In  summer  then  in  winter  (more  water  vapor).  Howe^'er,  the 
slop#  of  the  empirical  line  la  a little  too  steep  for  the  line  to 
go  through  the  origin.  The  derivation  of  Eq.  (23)  Included  some 
approxlmatiuns,  which  must  be  studied  lurcher. 

For  correcting  range  data,  it  is  not  actually  the  bending 
sngls  that  is  nsodad  but  tliC  effect  of  the  signal  path  curva.ure 
on  maasursd  range.  Figure  6,  obtained  from  ray  tracing,  shows  this 
affect  lor  signals  arriving  at  the  same  low  elevation  angle  (3®) 
under  various  wsochcr  conditions.  The  data  sample  includes  data 
for  the  first  half  of  January  and  July,  respectively  (same  sample 
ns  in  Fig.  7),  The  elf  art  of  bending  on  measured  range  is  plotted 
against  surface  refract! vity.  There  is  a fairly  linear  relation 
hut  with  connldernbla  scatter,  doubtless  due  to  refractlvlty  pro- 
file Irregular Jtlen.  Figures  7 and  8 show  obvious  correlation,  but 
N^yrfacf  ® better  predictor  of  the  bending  angle  than  of  fro 
range  effect.  It  Is  possible  that  Fgurfacc  ^ better  pre- 
dictor of  these  effects  than  fot  the  dry  component  alone 
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Fig,  7 Total  Angular  Bending  of  Signal  Path  versus  Surface  Refractivity 
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Fig.  8 Range  Effect  of  Signal  Path  Bending  versus  Surface  Refractivity 
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(as  it  is  for  high-angle  effects) , but  this  is  not  likely  to  be 
true  for  the  radio  case.  The  matter  is  being  studied  further. 

Equation  (23)  provides  a means  of  predicting  the  bending 
angle  from  surface  data  at  any  angle,  if  a value  is  postulated 
for  the  ratio  (R)  of  Eqs.  (19)  and  (20).  Figures  9 and  10  (opti- 
cal and  radio,  respectively),  show  this  predicted  bending  as  a 
function  of  elevation  angle  (E  = (ir/2)  - 8)  for  an  observed  atmo- 
sphere, and  also  the  bending  determined  by  ray  tracing  through 
the  same  atmosphere.  Three  different  values  of  R are  used  for 
theoretical  prediction  in  each  case.  The  differences  among  them 
are  insignificant  at  high  elevation  angles  but  can  be  seen  below 
10* . In  each  case,  the  ray-traced  values  of  x follow  one  of  the 
theoretical  predictions  down  to  approximately  3®  elevation  and 
depart  from  theory  below  that.  However,  different  values  of  R 
are  needed  in  the  optical  case  than  in  the  radio  case.  The 
larger  value  that  is  needed  in  the  radio  case  implies  a lower 
"center  of  gravity"  for  the  radio  bending  effect  than  for  the  op- 
tical bending  effect;  this  would  be  expected  because  of  the  water 
vapor  distribution. 

It  may  be  possible  to  develop  an  empirical  weighting  func- 
tion for  R so  as  to  make  a prediction  valid  at  still  lower  an- 
gles. Climatic  and  weather  effects  on  the  best  value  of  R also 
need  to  be  investigated. 

We  have  no  theoretical  expression  for  the  effect  of  path 
curvature  on  range  error  estimation  in  the  real  atmosphere  as  a 
function  of  bending  angle.  Figure  11  shows  the  effect  computed 
by  ray  tracing.  If  the  entire  path  were  an  arc  of  a single  cir- 
cle, the  curvature  effect  on  range  would  vary  as  x^,  but  this  is 
not  the  case.  Figure  11  shows  the  path  curvature  effect  as  a 
function  of  x for  one  observed  atmosphere.  The  observed  line  on 
the  log-log  paper  is  almost  straight  and  its  slope  indicates  that 
the  curvature  effect  in  that  atmosphere  was  proportional  to  x^ 
where  2 < K < 3 and  K probably  has  a value  not  far  from  2.5. 

No  statistical  study  of  this  relation  has  been  made  as  yet.  Fur- 
ther work  is  planned . 


MODEL  PREDICTIONS  AT  LOW  ANGLES 

In  developing  an  oblique-angle  correction  from  the  tropo- 
sphere model,  it  was  assumed  that  bending  of  the  signal  can  be 
neglected.  Range  effects  obtained  from  the  model  should  there- 
fore be  compared  to  values  obtained  by  numerical  integration  along 
a straight  line,  as  in  the  preceding  section.  If  they  can  be  made 
to  match  at  low  as  well  as  at  high  angles,  it  will  still  be  nec- 
essary to  apply  a curvature  correction. 
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Fig,  10  Radio  Signal  Path  Bending,  Theoretical  and  Ray  Traced 
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The  possibility  should  also  be  considered  that  the  model  or 
Its  parameters  could  be  adjusted  to  match  results  computed  along 
the  curved  path  at  low  and  at  high  angles.  In  this  case,  no  sepa- 
rate curvature  correction  would  be  needed,  a result  that  is  pref- 
erable but  less  likely  to  be  achieved.  We  first  must  learn  how 
actual  model  predictions  compare  with  straight-line  and  curved- 
path  numerical  Integration  results. 

Figures  12  and  13  show  the  difference  between  predictions 
from  the  model  using  surface  data  and  numerical  Integration  along 
a straight-line  path  using  balloon  data.  The  figures  show  a 
winter  and  a summer  sample.  Two  values  of  the  water-vapor  height 
parameter  for  the  model  are  used  for  each  month. 

The  discrepancy  between  the  model  and  the  observed  atmo- 
sphere varies  with  elevation  angle.  For  most  cases,  it  in- 
creases toward  the  horizon.  It  is  possible  to  match  the  two  ef- 
fects perfectly  at  the  zenith  for  any  reasonable  model  profile  by 
choice  of  parameters,  but  only  a perfect  match  of  refractivity 
profiles  (model  versus  observed)  would  match  effects  perfectly  at 
lower  angles  also.  Work  will  continue  on  improving  the  matching. 

The  lower-angle  error  that  is  introduced  by  the  present 
model  is  sensitive  to  the  choice  of  height  parameters,  as  can  be 
seen  in  the  figures.  The  parameter  for  dry-component  height,  in 
the  latest  form  of  the  model  (used  in  these  figures),  is  a linear 
function  of  surface  temperature.  No  weather  adjustment  is  now 
used  for  the  wet-component  height,  but  low-angle  data  may  provide 
criteria  for  making  such  an  improvement.  Further  work  along  this 
line  is  planned. 

Comparison  of  Figs.  12  and  13  with  Fig.  3 shows  that  sig- 
nal path  bending  is  a larger  contributor  to  prediction  error  at 
low  angles  than  is  model  inaccuracy.  The  total  prediction  error 
includes  both  of  these.  At  some  angles  and  with  some  parameter 
values  the  two  range  effects  may  have  opposite  signs,  so  that  a 
tradeoff  might  be  possible.  Figures  lA  and  15  show  winter  and 
summer  samples,  respectively,  of  the  difference  between  the  model 
prediction  and  the  actual  tropospheric  effect  computed  along  the 
curving  signal  path.  The  values  shown  are  one-month  means,  and 
error  bars  show  the  standard  deviation  from  the  mean.  These  are 
the  prediction  errors  that  should  be  expected  after  the  present 
model  correction  is  made.  The  mean  errors  for  winter  and  summer 
are  not  too  different,  but  the  scatter  is  much  worse  in  summer 
because  of  the  large  and  variable  water-vapor  content. 
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EFFECT  OF  TROPOSPHERIC  CORRECTION  ON  STATION  POSITION 

When  satellite  doppler  data  are  used  to  determine  the  posi- 
tion of  a tracking  station,  the  uncorrected  tropospheric  effect 
introduces  an  error  into  the  computed  station  position.  The  ap- 
parent range  to  the  satellite  is  too  great  at  every  point  of  the 
pass,  but  this  error  is  very  much  greater  at  the  horizon  than  at 
the  point  of  closest  approach.  The  time  derivative  of  the  ob- 
served doppler  shift  is  therefore  too  great,  and  the  station  ap- 
pears to  be  closer  to  the  orbit  than  it  really  is.  The  effect  is 
greater  for  low-angle  than  for  high-angle  satellite  passes.  Only 
the  apparent  range  is  affected  if  tracking  data  are  symmetrical 
about  the  point  of  closest  approach;  if  not,  the  apparent  station 
position  may  also  be  shifted  parallel  to  the  orbit. 

Figure  16  shows  the  range  effect  only.  The  range  error, 
with  no  tropospheric  correction,  determined  for  each  of  a set  of 
passes  at  the  Maine  OPNET  station  is  plotted  as  a function  of 
elevation  angle  of  the  satellite  at  closest  approach.  Each  point 
represents  a satellite  pass.  The  entire  computation  was  then  re- 
peated with  the  tropospheric  correction  included.  Both  sets  of 
results  are  shown  in  Fig.  16.  No  weather  data  were  recorded,  and 
statistical  values  were  used  as  a basis  for  the  computed  correc- 
tion. 


Range  errors  without  tropospheric  correction  varied  from 
approximately  12  m for  high-angle  passes  to  55  m for  some  lower- 
angle  passes.  Use  of  the  correction  left  little  range  error  in 
the  high-angle  passes,  but  the  low  passes  were  apparently  over- 
corrected. This  might  be  related  to  height  parameters  (cf.  Figs. 
12  and  13),  to  the  neglected  bending  effect  (Fig.  3),  or  to  the 
use  of  statistical  weather  data  instead  of  observed. 

Figure  17  shows  the  same  kind  of  data  but  for  a different 
location,  the  Hawaii  OPNET  station.  Again,  the  high-angle  passes 
were  corrected  adequately.  The  correction  for  the  lowest  pass, 
however,  was  too  small,  in  contrast  to  Fig.  16.  Again,  there  is 
a question  about  height  parameters  and  surface  weather  data. 

Figures  16  and  17  include  data  down  to  an  elevation  angle 
of  1°  when  available.  The  navigation  procedure  was  repeated  with 
a data  cutoff  angle  of  5°,  with  no  appreciable  change  in  the  re- 
sults. The  reason  is  not  certain,  but  it  is  possible  that  most 
of  the  passes  were  actually  not  tracked  much  lower  than  5°.  It 
would  be  desirable  to  continue  this  study  in  more  detail. 

The  lack  of  surface  meteorological  data  for  the  computa- 
tions that  produced  Figs.  16  and  17  means  that  the  correction  was 
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Fig.  16  Effect  of  Tropospheric  Refraction  on  Apparent  Tracking  Station  Position, 
Winter  Harbor,  Maine 


Effec\  ol  Tropospheric  Refraction  on  Apparent  Tracking  Station  Position, 
Wahiawa,  Hawaii 
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based  on  seasonal  values  of  the  refractlvlty , not  on  actual  values 
at  the  time.  This  lack  also  precludes  the  use  of  the  latest  formu- 
lation for  h(]  (Refs.  3 and  4).  The  use  of  observed  surface  data 
could  not  be  expected  to  change  the  navigation  results  by  more 
Chan  a few  meters,  even  for  low  passes.  However,  the  Importance 
of  a few  meters  has  Increased  In  recent  years. 

The  positions  of  Figs.  16  and  17  were  recalculated  using  a 
larger  value  for  h„  than  before,  12  km  instead  of  10  km.  The  cor- 
rections were  a little  larger  chan  before,  as  expected,  and  shifted 
apparent  station  positions  in  the  slant  range  coordinate  by  a few 
centimeters  for  the  very  high  passes  and  about  1 m for  the  low 
passes.  Changes  in  the  dry  component  would  have  more  effect,  but 
no  further  experimentation  has  beer,  carried  out. 
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4.  SUMMARY  OF  FRESENT  STATUS  AND  FUTURE  PLANS 


The  purpose  of  this  effort  is  to  make  poealblc  the  use  of 
satellite  range  and/or  range-rate  data  for  precise  studies  at 
lower  angles  than  are  now  practicable.  To  this  end,  the  tropo- 
spheric refraction  correction  must  be  as  accurate  as  possible. 

A correction  that  was  based  on  high-elevation-angle  con- 
siderations was  the  result  of  earlier  studies  (Refs.  2 through  3). 
Although  that  correction  can  be  used  at  any  angle,  its  accuracy 
deteriorates,  by  a moderate  percent,  toward  the  horizon.  This 
report  describee  the  progress  made  to  date  in  learning  more  about 
the  low-angle  effects  and  how  to  Include  them  in  an  improved  tro- 
pospheric correction.  Rending  of  the  signal  path,  neglected 
earlier,  is  now  taken  into  account.  Also,  adjustment  of  the 
mathematically  modeled  correction  for  better  correspondence  with 
the  real  atmospliere  is  under  study. 

Plans  for  further  study  may  be  summarized  under  the  follow- 
ing broad  categories.  Some  of  the  Itsmo  hevs  been  mentioned  in 
the  preceding  discussion. 


UPDATING  OF  THE  DATA  BASE 

Most  of  the  meCeorologicsl  dsts  on  hand  are  for  the  yssr 
1967,  the  most  recent  data  available  when  the  detelled  etudy  of 
balloon  data  was  begun.  Only  smell  amounte  of  older  and  newer 
data  have  been  obtained.  Balloon  date-gathering  technlquoe  hovi* 
Improved  in  the  intervening  years.  In  particular,  there  are  now 
better  Instruments  for  measuring  relative  humldlry.  Allltough  the 
water-vapor  component  accounts  for  only  a smoll  frertlon  of  atmo- 
spheric refractlvity,  it  is  Important  because  It  in  t utiivnaih}  w for 
most  of  the  variation  and  the  difficulty  in  sstlmating  the  effert 
from  surface  data.  Study  of  more  roceiil  balloon  dels  la  Indigstad, 

There  is  a further  possibility  of  a aids  benefit  heie,  Kk- 
tension  of  the  d/ita  base  in  time  might,  alao  provide  IMiiin,  at  Ivaai, 
of  a solar-cycle  effect  on  the  atmoaplierc , (It  la  also  of  Iniar- 
est  thot  the  cost  of  obtaining  balloon  data  from  the  Nat  loii/il  Cli- 
matic Center  hns  decreased  hfraua**  of  lni|irnved  proceduroa  at  l he 
Center. ) 
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MORE  RAY-TRACING  STUDIES 

The  present  troposphere  model  should  be  compared  with  addi- 
tional observed  data  covering  different  geographical  locations. 

Jhe  effects  of  varying  the  model  parameters  should  be  further  in- 
vestigated. The  present  parameters  are  designed  for  high  eleva- 
tion angles.  A method  should  be  sought  for  adjusting  the  model  to 
fit  propagation  at  all  angles,  low  as  well  as  high.  Possibly  a 
multiplicative  correction  factor  can  be  found,  e.g.,  a function  of 
the  elevation  angle  (F(E)),  a "correction  for  the  correction," 

This  factor  might  provide  the  correction  for  agreement  with 
hypothetical  straight-line  propagation  through  an  observed  atmo- 
sphere (in  which  case,  a correction  for  the  curvature  effect  must 
be  provided  separately),  or  it  might  include  the  curvature  effect. 
The  choice  itself  is  a matter  for  study.  Either  case  calls  for 
more  study  of  the  bending  effect  with  regard  to  climatic  and 
weather  variations. 


MORE  NAVIGATION  STUDIES 

Navigation  results  should  be  examined  for  a larger  sample 
of  stations  and  satellite  passes.  Preferably,  locally  observed 
surface  weather  data  should  be  provided  with  the  tracking  data  and 
used  with  the  model.  The  effects  of  varying  the  model  parameters 
and  of  changing  the  data  cutoff  angle  should  be  studied  and  com- 
pared with  expectations  based  on  ray-tracing  results. 

If  a new  correction  factor  (F(E))  for  improving  the  tropo- 
spheric correction  is  developed  by  ray-tracing  techniques,  the 
final  evaluation  should  make  use  of  satellite  data. 

Even  if  local  balloon  data  were  available  for  the  time  of  a 
satellite  pass,  predictions  based  on  ray  tracing  could  not  be  ex- 
pected to  agree  precisely  with  real  atmospheric  effects  on  the 
satellite  signal.  A balloon  gives  a profile  of  the  refractivity 
that  we  have  assumed  to  be  vertical,  but  the  ascent  is  in  fact  not 
vertical.  Alsu,  we  are  not  able  to  take  into  account  irregulari- 
ties such  as  clouds  and  updrafts.  The  idea  of  a "layered"  atmo- 
sphere is  an  approximation,  while  the  satellite  signal  traverses 
the  real  atmosphere.  Thus  use  of  the  model  with  navigation  data 
provides  the  acid  test.  As  long  as  use  of  the  model  results  in 

any  systematic  navigation  errors,  further  improvement  should  be 
possible. 


im  /OHN»  MO»«'N»UNrVl«li'* 

A^UCO  PHYflCS  LABOSATOflV 


ACKNOWLEDGMENT 


Th«  cnll*bor*tlon  of  Harry  K.  Utterback  In  developing  the 
ray-Craclng  program  and  making  the  numerical  studies  baaed  on  It 
la  acknowledged  with  appreciation.  The  navlfstlon  program  that 
was  uaed  In  the  Later  section  of  the  work  was  developed  by  the 
Analyala  and  Compucstlon  Croup  of  the  Space  Development  Depart- 
nanc  of  this  Laboratory.  Thanks  are  due  tc  L.  Lee  Pryor  for  mak- 
ing epeclal  modifications  In  the  program  and  handling  the  compu- 
tations Involved. 


39 


ThJ  JOHNS  HOPKINS  UNIVenSITV 

APPLIED  PHYSICS  LABORATORY 
laurcl.  marylano 


REFERENCES 


1.  E.  K.  Smith,  Jr.,  and  S,  Weintraub,  "The  Constants  in  the 
Equation  for  Atmospheric  Refractive  Index  at  Radio  Frequen- 
cies," Proc.  I.R.E.,  Vol.  41,  No,  8,  August  1953,  pp.  1035- 
1037. 

2.  H.  S.  Hopfleld,  "Two-Quartlc  Tropospheric  Refractlvity  Pro- 
file for  Correcting  Satellite  Data,"  J.  Geophys.  Res. . 

Vol.  74,  No.  18,  20  August  1959,  pp.  4487-4499. 

3.  H.  S.  Hopfleld,  "Tropospheric  Effect  on  Electromagnetically 

Measured  Range:  Prediction  from  Surface  Weather  Data," 

Radio  Scl. . Vol,  6,  No.  3,  March  1971,  pp.  357-367. 

4.  H.  S.  Hopfleld,  "Tropospheric  Range  Error  Parameters;  Fur- 
the.  Studies,"  APL/JHU  CP  015,  June  1972  (also  published  as 
Goddard  Space  Flight  Center  Preprint  X-551-72-285,  August 
1972). 

5.  S.  M.  Ylonoulls,  ’’Algorithm  to  Compute  Tropospheric  Refrac- 
tion Effects  on  Range  Measurements,"  J.  Geophys.  Res..  Vol. 
75,  No.  36,  20  December  1970,  pp.  7636-7637. 

6.  G.  Joos,  Theoretical  Physics,  translated  by  Ira  M.  Freeman, 
Hafner  Publishing  Co.,  New  York,  1934. 

7.  J.  E.  Freehafer,  "Geometrical  Optics,"  Vol.  13,  Chap.  2, 
Propagation  of  Short  Radio  Waves,  D,  E,  Kerr  (Ed.),  M.  I.  T. 
Radiation  Laboratory  Series,  McGraw-Hill  Book  Co.,  1951,  and 
Dover  Publications,  New  York,  1965. 

8.  W.  M.  Smart,  Text-Book  on  Spherical  Astronomy,  Cambridge 
University  Press,  Cambridge,  England,  1962. 

9.  B.  R.  Bean  and  F.  J.  Dutton,  "Radio  Meteorology,"  National 
Bureau  of  Standards  Monograph  92,  1 March  1966, 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY* 

Th«  work  Ttportad  In  TG  1291  was  done  under  Navy  Contract  N00017-72-C-4401.  This  work  is  related  to  Task  TIC, 
which  la  supported  by  Defense  Mapping  Agency  Topographic  Center,  52400. 


ORGANIZATION 


DEPARTMENT  OF  DEFENSE 

enae  Napping  Agency  Headquarters 
ante  Napping  Agancy  Topographic  Ctr. 


Deparfent  of  the  Navy 


LOCATION 


ATTENTION 


ATAIRSTSCOH 

AVSEASYSCOK 

aval  Oceanographic  Office 
laval  Research  Laboratory 


Naval  Surface  Weapons  Ctr. 

NAVPRO 

Deoartnant  of  the  Air  Force 

AF  Western  Test  Range 
AFCRL,  Hanscoo  Field 

1160th  Support  Group, ' Bolling  AFB 

Depar tment  of  the  Arm 


A.C.I.C.  Detachment  No.  1 

GOVERNMENT  AGENCIES 
NAS\  Goddard  Space  Flight  Center 


MASA  Wallops  Station 
National  Bureau  of  Standards 
NOAA  Geodetic  Research  and  Development 
Laboratories 

Smithsonian  Astrophyslcal  Observatory 


SCHOOLS 

New  Mexico  State  University  Physical 
Science  Laboratory 

Massachusetts  Institute  of  Technology 
University  of  Maryland,  Department  of 
Physics  and  Astronomy 


Ohio  State  University,  Department  of 
Geodetic  Science 
University  of  Texas 


Washington,  DC 
Washington,  DC 


Alexandria,  VA 


Washington,  DC 
h’ashlngton,  DC 
Washington,  DC 
Washington,  DC 


Dahlgren,  VA 
Silver  Spring,  MD 


Vandenberg,  CA 
Bedford,  MA 

Washington,  DC 


Arlington,  Va. 


Greenbelt,  MD 


Wallops  Island,  VA 
Boulder , CO 

Rockville,  MD 
Cambridge,  MA 


University  Park,  NM 
Cambridge,  MA 

College  Park,  MD 


Columbus,  OH 
Austin,  TX 


Charles  Martin 

41300 

52200 

52220 

52400 


Library,  AIR-50174 
Library,  SEA-09G3 
D.  Des  jardins,  3541 
J.  Bulsson,  5165 
L.  D.  Breetz 
L.  Choy,  7112  CH 
R.  Eislnger,  7918 
DK-10 


C.  Buchelt 
R.  W.  Lee 
E.  Altschuler 
CAPT  P.  Golluccl 


J.  W.  Marini 

P.  E.  Schmid 

J.  Herbert 

D.  M.  LeVine 

S.  Gubin 

C.  D.  Leltao 

M.  C.  Thompson,  Jr 

B.  Witte 

G.  C.  Welffenbach 

C . G . Lehr 

M.  Gaposchkin 

M.  R.  Pearlman 

N.  C.  Mathur 


D.  Martin 
I.  I.  Shapiro 

C.  0.  Alley 

S.  K.  Poultney 

D.  C.  Currie 
I.  I.  Mueller 

F.  A.  Fajamlrokun 
A.  Tucker 


'Inillsl  dUtKbution  of  Ih»  document  within  the  Applied  PhyfiM  Leboretory  h»i  been  made  in  eccordanca  with  a lljt  on  file  in  the  APL  Technical  Publications  Group. 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY* 

TG  1291 
-2- 


ORGANIZATION 

LOCATION 

ATTENTION 

No.  of 
Copies 

CONTRACTORS 

California  Institute  of  Technology,  Jet 

G.  C.  C.  Chao,  CPB-lOii 

1 

Tropulslon  Uiboratory 

Pasadena,  CA 

P.  F.  MacDoran 

1 

V.  J.  Ondraslk 

1 

R.  Jaffe 

1 

T.  Otoshi 

1 

Joint  Institute  for  Laboratory  Astro- 

physics,  University  of  Colorado 

Boulder,  CO 

P.  L.  Bender 

1 

Lincoln  Laboratory,  Massachusetts 

Institute  of  Technology,  Millstone  Hill 

Radar 

Lexington,  MA 

A.  Freed 

1 

Analytical  Sciences  Corp. 

Reading,  MA 

W.  Larlmore 

1 

DBA  Systems 

Lanham,  MD 

G.  Klebuzlnskl 

1 

Mitre  Corp. 

McLean,  VA 

B.  D.  Elrod 

1 

General  Electric  Co. 

Utica,  NY 

T.  Godbey,  MD  309 

1 

Interstate  Electronics  Corp. 

Anaheim,  CA 

Tech.  Library,  Prolantonlo 

1 

Offshore  Navigation  Inc. 

Harahan , LA 

E.  H.  Christy 

1 

Satellite  Positioning  Corp. 

Enc  ino , CA 

D.  Ozdes 

1 

Wolf  Research  and  Development 

Rlverdale,  MD 

1 

Magnavox  Corp. 

Torrance,  CA 

T.  Stansell 

1 

Aerospace  Corp. 

El  Segundo,  CA 

H.  M.  Wachowski 

1 

Los  Angeles,  CA 

S.  Leonard 

1 

Collins  Radio  Co. 

Dallas,  TX 

J.  R.  Hoelzel 

1 

Boeing  Scientific  Research  Lab. 

Seattle,  WA 

W.  G.  Tank 

1 

Honeywell  Inc.,  Gov't  and  Aeronautical 

Products  Division 

Meanneapolls,  MN 

R.  J.  Kell,  Mail  Sta.  A3607 

1 

Systems  Group  of  TRW  Inc.,  Washington 

Operations 

Washington,  DC 

W.  H,  Lewellen 

1 

Rand  Corp. 

Santa  Monica,  CA 

J.  L.  Bower 

1 

RCA 

Morristown,  NJ 

R.  Mitchell 

1 

FOREIGN 

Geodetic  Survey  of  Canada 

Ottawa,  Canada 

J.  D.  Boal 

1 

J . Kouba 

1 

H.  E.  Jones 

1 

Atlantic  Oceanographic  Lab., 

Dartmouth,  Nova  Scotia, 

Bedford  Institute 

Canada 

D.  Wells 

1 

York  University,  Faculty  oi'  Science 

Downsvlew,  Ontario,  Canada 

R.  B.  Langley 

1 

University  of  New  Brunswick,  Department 

Fredericton,  New  Brunswick, 

of  Surveying  Engineering 

Canada 

M,  M,  Nassar 

1 

National  Research  Council  of  Canada, 

Division  of  Physics 

Ottawa,  Canada 

J.  Saastamolnen 

1 

Canadian  General  Electric  Co.  Research 

Dept. 

Toronto,  Ontario,  Canada 

M.  R.  Holmes 

1 

Brigadier  G.  Bomlord 

Sutton  Courtenay,  Oxon,  England 

1 

Centre  National  d 'Etudes  des 

Tll^communlcatlons 

Issy-les-Moulineaux,  France 

A.  E.  Azoulay 

1 

Centre  National  d' Etudes  Spatiales 

Brdtigny  sur  Orge,  France 

F.  Mouel 

1 

M.  P.  Guerlt 

1 

Observatolre  Royal  du  Belgique 

Brussels,  Belgium 

P.  Paquet 

1 

Instltut  fUr  Angewandte  Ceodbsle 

Frankfurt,  Germany 

P.  Wilson 

1 

K.  Nottarp 

1 

Astronomical  Institute  of  the 

Czechoslovak  Academy  of  Sciences 

Ondrejov,  Czechoslovakia 

M.  Prokeil 

1 

Meteorological  Service,  Division  for 

Research 

Bet  Dagen,  Israel 

A.  Kanes 

1 

National  Physical  Laboratory.  Radio 

Propagation  Unit 

New  Delhi,  India 

Y,  V.  Semayajulu 

1 

National  Physical  Laboratory,  Radio 

Science  Division 

New  Delhi,  India 

S.  B.  S.  S.  Sarma 

1 

Ministry  of  Defense,  Hlm.slayan  Radio 

Propagation  Unit 

Raipur,  Dchra  Dun,  India 

E.  Bhagir.Ttha  Rao 

1 

Space  Applications  Center,  Geodesy 

Jodhpur  Tekra,  Ahmedabad, 

Sec 1 1 on  (RSA) 

Ind  in 

M.  R.  Sivaraman 

1 

National  M;.ipptn^' 

C.snbcrr.s  City,  Australi.s 

B.  P.  L.nmbcrt 

1 

•Initial  di»Inbu»ion  of  thi5  document  within  the  Applied  Physics  Laboratory  hes  been  made  in  accordance  with  a list  on  file  in  the  APL  Technical  Publications  Group. 


INITIAL  DISTRIBUTION  EXTERNAL  TO  THE  APPLIED  PHYSICS  LABORATORY^ 


TG  1291 
-3- 


ORGANIZATION 


LOCATION 


ATTENTION 


FOREIGN  (cont'd) 

University  of  New  South  Wales, 
School  of  Surveying 
Weapons  Research  Establishment, 
Space  Research  Group 
Australian  Army  Staff 


Kensington,  N.S.W. , Australia 

Adelaide,  S.A. , Australia 
Washington,  DC 


A.  Stolz 
I.  D.  Lloyd 

K.  J.  W.  Lynn 
LT  D.  McLuskey 


Initnl  diitrlbution  of  Ihu  doruowol  within  the  AppluKf  Phyitct  Lubomtofy  h»»  h*«n  mad*  In  accordimca  with  a litt  on  fils  In  the  APL  Technical  Piihlicjitiont  Group. 


